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1.0 SUMMARY 


The program L218 (LOADS) is structured as five overlays, one main overlay and four 
primary overlays. Input into the program is made via cards and magnetic files (tapes or 
disks). Output from the program consists of printed results and magnetic files 
containing load coefficient matrices for use in either L219 (EQMOD) or L221 (TEV156). 

Although L218 (LOADS) serves as a module of the DYLOFLEX system, it can be 
operated as a stand alone program. Subroutines used by L218 include routines 
embedded in the program code and routines obtained from the standard FORTRAN and 
the DYLOFLEX libraries. 



2.0 INTRODUCTION 


The computer program L218 (LOADS) was developed for use as either a standalone 
program or as a module of a program system called DYLOFLEX (see fig. 1) developed 
for NASA under contract NAS1-13918 (ref. 1). Because of the DYLOFLEX contract 
requirements defined in reference 2, a program was needed to calculate dynamic load 
coefficient matrices for use as sensors in active control analyses and/or for use in 
calculating design loads. 

The objective of this volume is to aid those persons who will maintain and/or modify the 
program in the future. To meet this objective, the following items are defined: 

• Program design and structure 

• Overlay purpose and description 

• Input, output and internal data base descriptions 

• Extent of checkout 
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Figure 1. — Dyiof/ex Flow Chart 








3.0 PROGRAM DESIGN AND STRUCTURE 


The program is structured as a system of five overlays (fig. 2): 


Main overlay (L218,0,0) 
Primary overlay (L2 18,1,0) 
Primary overlay (L218,2,0) 
Primary overlay (L218,3,0) 
Primary overlay (L218,4,0) 


L218vc 

RGEN 

AVD 

NPLDS/PLDS 

VBMT 


The main overlay L218vc controls reading of general data cards by the primary overlay 
(1,0) RGEN, and calls the proper primary overlay to perform the execution requested. 
The main overlay does not read input cards. 


The (1,0) primary overlay RGEN reads the general input data and the module cards 
($AVD, $NPLDS, $PLDS, or $VBMT), which selects the primary overlay for execution. 


The (2,0) primary overlay AVD reads the AVD card and tape input data, performs the 
AVD calculations, and writes the AVDTAP output tape. 

The (3,0) primary overlay NPLDS/PLDS reads the NPLDS (or PLDS) card and tape 
input data, performs the NPLDS (or PLDS) calculations, and writes the NPTAP (or 
PTAP) output tape. 

The (4,0) primary overlay VBMT reads the VBMT card and tape input data, performs 
the VBMT calculation, and writes the LTAP output tape. 

Although L218 serves as a module of the DYLOFLEX system, it can be operated as a 
standalone program. When the program is run by itself, it becomes the user’s 
responsibility to generate input data in the format required by L218; see figure 3 for the 
external file requirements and volume I of this document for the file contents and 
formats. 


This program requires subroutines that are not part of the L218 code. Some routines are 
automatically obtained from the standard FORTRAN library when the program is 
loaded. Others, however, are stored in the DYLOFLEX alternate subroutine library, 
which must be declared at the time of loading. Subsequent sections describe each 
overlay separately and contain tables displaying the routines called and the library in 
which they are located. All subroutines are single entry point routines. 

This volume describes the program in a macro sense. A more detailed discussion 
appears in the comments contained in the program source code. Each routine contains a 
preface describing the routine’s purpose, author, analytical steps, modification history, 
input data, output data, glossary of variables, and list of other routines called. 
Embedded within the executable code are comments labeling each section and 
explaining logical branches. 
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Figure 2. — Overlay Structure of L2 18 (LOADS) 










Figure 3. — Input 'Output of L21 8 Overlays 












3.1 OVERLAY (218,0,0) - L218vc 


The main overlay of L218 is itself named L218vc where v is a letter indicating the 
program version and c is an integer number indicating the correction number that 
applies to the v version. 

Purpose of L218vc 

Overlay L2l8vc directs the execution of the primary overlays and aids communication 
between primary overlays via labeled common blocks. 

Analytical Steps of L218vc 
Overlay L2l8vc performs its task in the following steps: 

1. Overlay (L218,l,0), RGEN, is called to read general input and/or a keyword card 
$AVD, $NPLDS, $PLDS, or $VBMT. 

2. L218 determines the next overlay to call. 

a. If the keyword is $AVD, jump to step 3. 

b. If the keyword is $NPLDS, jump to step 4. 

c. If the keyword is SPLDS, jump to step 4. 

d. If the keyword is $VBMT, jump to step 5. 

e. If the keyword is $QUIT, jump to step 6. 

3. Overlay (L218,2,0) is called to read input data and perform calculations for AVD 
(acceleration, velocity, displacement). When finished, program control returns to 
step 1. 

4. Overlay (L218,3,0) is called to read input data and perform calculations for 
NPLDS/PLDS (net panel loads or panel loads). When finished, program control 
returns to step 1. 

5. Overlay (L218,4,0) is called to read input data and perform calculations for VBMT 
(shears, bending moments, and torsions). When finished, program control returns 
to step 1. 

6. L218 program stop. 

The macro flow chart of this overlay is shown in figure 4. The subroutines called are 
displayed in table 1. 

I/O Devices of L218vc 

There is no I/O performed in the (0,0) overlay. It is all performed in the primary overlays. 
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Table 1. — Routines Called by L218vc 


OVERLAY 

(L2 1 8 # 0 r 

0) 

PRCGFAM 

L2 1 3 vc 


RGEN 

OVERLAY 

(L218, 1 ,0) 

AVD 

OVERLAY 

(L2 18,2,0) 

NPLDS 

OVERLAY 

(L2 18,3,0) 

VBMT 

OVERLAY 

(L218,4,0) 
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3.2 PRIMARY OVERLAY (L218,l,0) - RGEN 
Purpose of RGEN 

The L218 (LOADS) program’s first primary overlay is named RGEN. RGEN reads 
general input data from cards and may write a J-MATRIX and a MASS MATRIX tape. 
Upon reading a keyword defining the LOADS module to be executed, control is returned 
to L218vc. 


Analytical Steps of RGEN 

RGEN performs its task in the following steps: 

1. If this is not the first call to RGEN, jump to step 4. 

2. The header matrix array (IQLTAP) is initialized. The subroutine PRGBEG is 
called to place the program header on the printed output. 

3. A data card is read. It must begin with $LOAD to assure that the card input file is 
correctly positioned. If it does not contain $LOAD the fatal error counter is 
incremented and the program tries to process additional cards; however, module 
execution will not occur. 

4. A program directive card is read, printed, interpreted, and acted upon according to 
the following conditions: 

a. If the keyword is $TITLE, begin step 4 again. 

b. If the keyword is $GEN, jump to step 5. 

c. If the keyword is $AVD, $NPLDS, $PLDS, or $VBMT, jump to step 11. 

d. If the keyword is $END, jump to step 4. 

e. If the keyword is $QUIT, jump to step 11. 

5. A general input directive card is read, printed, interpreted, and acted upon as 
follows: 

a. If the keyword is JMAT, jump to step 6. 

b. If the keyword is MASS, jump to step 7. 

c. If the keyword is TRAN, jump to step 8. 

d. If the keyword is MAXSUR, jump to step 9. 

e. If the keyword is AVDTAP, EOMLOD, JTAPE, LTAPE, MASSTP, NPTAP, 
IPTAP, or SATAP, jump to step 10. 

f. If the keyword is $ — , jump to step 4a. 
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6. Subroutine RGEN2 is called. This subroutine reads the J-matrix data and writes it 
on JTAPE. Program control returns to step 5. 

7. Subroutine RGEN1 is called. This subroutine reads the MASS matrix data and 
writes it on MASSTP. Program control returns to step 5. 

8. The transformation order is set. Program control returns to step 5. 

9. The maximum number of surfaces is set. Program control returns to step 5. 

10. The appropriate tape name is changed. Program control returns to step 5. 

11. The appropriate keyword code is set. If the keyword was not $QUIT, jump to 
step 13. If it was a $QUIT card, jump to step 12. 

12. Copy the final general output from IOUT to IUTFIL. Call PRGEND to place the 
program termination message on the printed output. 

13. Return control to (L218,0,0). 

The intended order of card input for RGEN is: $LOAD card, $TITLE card, $GEN card, 

all GEN input for this execution run, $module card. 

The macro flow chart of this overlay is shown in figure 5. The subroutines called are 

displayed in table 2. 


I/O Devices of RGEN 

RGEN reads general loads card input and may write the J-matrix (JTAPE) and mass 
matrix (MASSTP) magnetic files. 

3.3 PRIMARY OVERLAY (L218,2,0) - AVD 
Purpose of AVD 

The L218 (LOADS) program’s second primary overlay is named AVD (acceleration, 
velocity, and displacement). AVD processes modal deflections (</>), modal slopes (</>^), and 
geometry data (BS, BBL, WL, 0 X , 0 y , 6 Z ) to generate the appropriate coefficient 
matrices required by L221 (TEV156) to calculate loads. The loads consist of 
translational and/or angular accelerations, velocities, and displacements at selected 
points on the structure. The axis system may be either reference or local (that system 
defined by the angular data ( 0 X , 0 y , 0 Z ) in the geometry input , or by 
(0 X #z) from card input. The resultant loads matrices are written on tape (AVDTAP). 


Analytical Steps of AVD 

AVD performs its task in the following steps: 

1. It initialize FETS for disk storage Files MERGE 1, MERGE 2, and MERGE 3. 

2. The subroutine RAVD is called to read the AVD card input. 

3. If matrices were card input, jump to step 9. 

4. The subroutine DISK is called to read geometry, modes, and SA array from 
(SATAP). 
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Figure 5. — Macro Flow Chart of Overlay ( L218 , 1,0) RGEN 
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Table 2. — Routines Called by RGEN 


OVERLAY (L218,1,C) 
PROGRAM RGEN 


FETACD 

+ 


FETDEL 

+ 


NAMFIL 

+ 


PPGBEG 

+ 


PRGEND 

+ 


RGEN1 


WRTETP + 

PGEN2 


WRTETP + 


+ indicates a routine in the DYLOFLEX alternate subroutine 
li fcrary 

All others are local to L2 18 (LOADS). 
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5. If interpolation is required, jump to step 8. 

6. If transformation (rotation) is not required, jump to step 9. 

7. The subroutine ROTATE is called to perform the required axis transformation 
(rotation); jump to step 9. 

8. The subroutine NTERP is called to perform the necessary interpolation. 

9. The subroutine MERGE is called to merge the matrices for this surface. 

10. If more surfaces are to be processed, jump to step 2. 

11. The subroutine AVDTAP is called to write the final tape (AVDTAP). 

12. The subroutine RETURNF is called to return all scratch files. 

13. Return control to (L2 18,0,0). 

The macro flow chart of this overlay is shown in figure 6. The subroutines called are 
displayed in table 3. 


I/O Devices of AVD 

AVD reads card input. Geometry, mode shapes, slopes, and SA arrays are read from 
SATAP, which is normally written by the DYLOFLEX interpolation program L215 
(INTERP), described in reference 3. 

Regular and general print options control the printed output. 

The coefficient matrices are written on the final output tape (AVDTAP) in a format 
acceptable to L219 (EQMOD) and L221 (TEV156) (refs. 4 and 5, respectively). 

3.3.1 AVD PROGRAMMING SPECIFICATIONS 

AVD Modal Data Interpolation Methods 

When the interpolation option is selected, axis transformation (rotation) is not allowed. 
Interpolation is performed in the local axis system to an arbitrary (x,y,z) reference axis 
coordinate. 

There are three interpolation methods available, one method for thin bodies and two 
methods for slender bodies. Their description is as follows. 

1. Interpolation of Thin Bodies (INTER=1): Where <t>\ ~ = = 0, and 

<t>Q y are obtained by calling AINTG, a subroutine in the DYLOFLEX library, with 
the appropriate SA array and multiplying the results by -1 to retain the proper 
sign convention of the modes. 


14 
























Table 3. — Routines Called by A VD 


OVERLAY 

(L218, 2,0) 


PROGRAM 

AVD 


AVDTAP i 

FETDEI + 



FETADD + 



WRTETP+ 


DISK i 

DISKI 

DISK5 READTP+ 


DISK2 

DISK 5 READTP + 


_ DISK5 

READTP+ 

FETADD+ 



MERGE 

MERGE 1 


NTEFP i 

NTEPP1 

AINTG+ 


NTERP5 



NTERP6 


RAVE 

RAVD1 « 

RAVD1A PAVD6 



RAVC6 


RAVD2 < 

RAVD6 



KOMSTR+ 



RAVD5B 


RAVD3 1 

RAVC5B 



RAVD6 


♦ Indicates a routine in the DYLOFLEX alternate subroutine library. 
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Table 3. — (Concluded) 


RAVD 


< RAVD4 

RAVD5 

RAVD5A 

RAVD6 

READTP* 


< RAVD5A 
RAVE5B 

< RAVD5A 
RAVE5E 


ROTATE 


ROT AT El 



2. Interpolation of Slender Bodies (INTER=2). This linear interpolation method 
between points is suitable for elastic axes that are straight and, preferably, 
parallel to the reference axis. To interpolate to a point x, given two structural node 
points I and 1+1 and their coordinate locations (BS, BBL, WL): 


I 


1+1 


A 

K. 


= (x - BS,)/(BS I+I - BSj) 

= A ( 0 yi+i'S) + 0 yi + LTT 0 

= A { 4 >z i+\~ <Pz i) + 4>z i + LT 


= A 


y x 


= A 


1 

-0e„ 

) + 00 v 

X I+1 

x i/ 

/ x l 


> 

- 00 

\ + 00 

yi+i 


I y\ 


-00 > 


z i+i 

1 ; 

1 z i 


where: 

LT = BBL X - BBLi 
LTT = WLx - WL r 

To interpolate to a point x, given only one structural node point I and its coordinate 
locations (BS, BBL, WL): 

I X 

• 


or 


I 


then: 


0v = 0 


X I 


<V X = 0 yi " LB 


- 0 Z j + LB 


Z I 


+ LTT 


LT 00, 
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where: 



LB = BS X - BSj 
LT = BBL X - BBLj 
LTT = WL X - WL, 


3. Interpolation of Slender Bodies (INTER=4). A linear interpolation method that 
interpolates to a node using rigid links attached to a reference node. To interpolate 
to a point x using only one structural reference node I: 


*y x ‘ S “ LB % + LTT % 

<p = 4> + LB 4>q + LT <pQ 

*-x L l yj xj 


0« = ^0, 


c i 


- 00 


y\ 


where I, LB, LTT, and LT are card input using the relationship defined in method 2. 
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Transformation of Axis 


To transform the modal data from the local axis to the inertia (or arbitrary) axis 
requires that each mode 0 X , 0 y , 0z»^0x» 0# v and0^ z at each node for the surface be 
transformed to the axis defined by 0 X , 9 y , 0 Z as follows: 


V 


— 




= 

T 



0>z 

and 









te y 

= 

1 



<Pe z 



^ z 

_ 






local axis 


local axis 


where T is the Euler transformation matrix for d x , 6 y , Q z (see vol. I, sec. 5 ). 

Note: All 0’s must be present; they are assumed zero if null. 

Selection of Modes 

In general, 0 X , 0 y , 0 Z , 0^ ,<t>n and 0a are available, but for a normal case only some 

x Vy v z 

will be selected. Keywords will be input to select the required modes and final Mj, M2, 
and M3 matrices. 

For the acceleration matrix (M3), either translation (0 X , 0 y , 0 Z ) and/or rotation ( 0 q , , 

<t>Q ,<f > 0 ) may be required. Furthermore, any combination of 0 X , 0 y , and 0 Z may be 

required (similarly for 0 ^ , 0 ^ and 0 ^ ). The same is also true for the M2 and Mj 
matrix. x y z 

As an example, if an M3 translational matrix for x, y, z is required, then M3 X , M3 y , 
M3 z will be selected. If the M3 rotational matrix for Z is required, then M3 Z (rotational) 
will be selected. 


20 


Merging of Matrices 

Matrices for a given surface are merged in the order x, y, z translational, x, y, z 
rotational. Any matrix that is null is omitted. If all of the matrices M3 were nonzero for 
surface IS, M3 for surface IS would merge as: 


M, 

J x 

1 

^3 V 


y 


Mo 

J z 



> 

Mo 



► 

y 


Mo 

-< 


translational 


rotational 


and similarly for M2 and Mj. This is further illustrated by the following examples: 

Assume the number of modes is three and that nodes 1 and 4 are selected from <f> Xt $ y , 
</>x and <f>Q , <f>0 for translational and rotational acceleration mode shapes, 

respectively. y 
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Then the merged matrix would appear as: 



If nodes 1 and 4 are selected from <f>y and <f> z for translational velocity only, M2 would 
merge as: 
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A VD’s Final Output Tape 

After all surfaces have been processed and the matrices merged, the final tape will 
contain the merged matrices [M^], and [M 3 ]. The maximum size of [Mj] is 

(100 x 70), similarly for [M 2 ] and [M 3 ]. 

To illustrate the format for final output, consider a load set of one surface with all 
matrices required (no zero matrices). Then: 
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These matrices would be written on the tape as: 


Header 

Mj 

m 2 

m 3 

EOF 


>Load set number 1 


rLoad set number 2 


EOF 


The physical size of the matrices M|, M2, and M3 are equal and are appropriately 
loaded with zero rows as follows. 

In general there may be several surfaces in a load set. Then: 


0 


0 


^surface 1 

0 


M ? 

surface 1 


0 

^ surface 1 


0 


0 

0 


0 


Mo 

■^surface 2 

0 

. [M 2 ] = 

m 2 

^surface 2 

[M 3 1 = 

0 

Mr , „ 

‘surface 2 


0 


0 

• 


. 


etc 

1 

CD 

r-+- 

O 

1 


etc 



where zero rows are added as required. 
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3.4 PRIMARY OVERLAY <L218,3,0) - NPLDS/PLDS 
Purpose of NPLDS/PLDS 

The L218 (LOADS) third primary overlay, NPLDS/PLDS, reads the specific card input 
data for net panel loads or panel loads. A net panel loads run will result in loads 
matrices written on tape NPTAP. A panel loads run will result in loads matrices 
written on tape PTAP. 

Analytical Steps of NPLDS/PLDS 

NPLDS/PLDS performs its task in the following steps: 

1. It initializes FETS for disk storage. 

2. The subroutine OPENMS is called to initialize a random access file, MERGB. 

3. The subroutine NPLDA is called to read card input data. 

4. The subroutine NPLDB is called to calculate M3 and/or read geometry for this 
surface (PLDS reads geometry but does not calculate M3). 

5. The subroutine NPLDD is called to calculate M4, M5, and 4 > for all frequencies. 

6. If more surfaces are to be processed, jump to step 3. 

7. The subroutine NPLDH is called to merge the matrices and write the final output 
tape NPTAP or PTAP. 

8. The overlay closes all FETS and returns scratch files. 

9. Return control to (L2 18,0,0). 

The macro flow chart of this overlay is shown in figure 7. The subroutines called are 
displayed in table 4. 


I/O Devices of NPLDS 

NPLDS reads card input. Geometry and modal data is read from tape (SATAP) as 
provided by L215 (INTERP) (ref. 3). Equations of motion data is read from tape 
(EOMLOD) as provided by L217 (EOM) (ref. 6). Mass matrix data is read from tape 
(MASSTP) as provided by (L218,l,0), RGEN. 

Note: Both SATAP and EOMLOD are required for either NPLDS or PLDS. 

Regular and general print options control the printed output. For net panel loads, the 
loads matrices are written on tape (NPTAP) in a format acceptable to L219 (EQMOD) 
(ref. 4) and L221 (TEV156) (ref. 5). 

For panel loads, the loads matrices are written on tape PTAP in a format acceptable to 
L219 (EQMOD) and L221 (TEV156). 
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Figure 7. — Macro Flow Chart of Overlay (L218,3,0) NPLDS/PLDS 














Table 4 . — Routines Called by NPLDS/PLDS 


OVERLAY (L218,3,0) 
PROGRAM NPLDS/PLDS 


FETADD+ 

FETDEL+ 

NPLDA 


NPLDB 


< FETADD+ 
NPLDA 1 

NPLDA2 

NPLDA3 

NPLDA6 
READTP+ 
/ NPLDB1 
NPLDB2 
PEADTP+ 


1 KOMSTR+ 
NPLDA6 
{ NPLDA 6 
READTP+ 
< NPLDA 4 
NPLDA6 


WRITMS* 


+ Indicares a routine in the DYLOFLEX alternate subroutine library. 
* Indicates a routine in the FORTRAN subroutine library. 
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Table 4. — (Concluded) 


NFLDD 


NPLDD1 

NPLDD2 

NPLDD4 


NPLCH 


OPENMS* 


•f FETDEL+ 
NPLDH1 
WR1ETP+ 


PEADTP+ 

/ NPLDD1 
NPLDD2 
NPLDD3 

WRITMS* 

- 

READMS* 


REACT P+ 

| AINTG+ 
PLATE* 
REA DTP* 


RETUPNF+ 


3.4.1 NPLDS/PLDS PROGRAMMING SPECIFICATIONS 


NPLDS (net panel loads) uses the modal deflection <f> z with the mass matrix to 
obtain the inertia forces (M3 matrix) on each subsurface. It also uses the aerodynamic 
force matrices generated in the Equations_of Motion program from which selected nodes 
are extracted to form the M 4 , M5, and $ matricies. The resulting matrices are in the 
local axis only. PLDS is identical to NPLDS except that the [M3] matrix is omitted. 
The following information is useful in understanding the program generation of the 
required matrices. 

For each requested surface and selected nodes: 

1. Calculate [M3] and write it on NPTAP. 

2. Form [M 4 ], [M5], and [</>] for all frequencies (k) and write them on NPTAP (also 
on PTAP if requested). 

To calculate M3: 

1. Read the mass matrix. 

2. Read scalars for [M3]. 

3. Read <f> z from disk (also geometry). 

4. Read nodes to be used. 

5. Calculate: 


m^ z 


1 1 


[M 3 ] = SCALR1 


m i^z 


il 


i^z; 


There is no contribution from </> x and <£ y . M3 is in the local axis only. 
To calculate [M 4 ], [M5], and [<£]: 

1. Read nodes. 

2. Read Fp]Jq), FpjJq), FpL(cEg) from the EOM tape. 

3. Read scalars for [M 4 ], [M5], and [<£]. 
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4. If the structural and aero node points are identical: 


[M 4 ] = SCALR2 


PL 


11 


[M 5 ] = SCALR2 


[0] = SCALR3 


where g = number of gust zones. 


"PL 


In 


l_ FPL il 

FpL inJ 

Fp l-I , 

FpL ln 

l_ FPL il 

FpL inJ 

Fp M i 

^lg 

Fpp „ 

Fp j 


CUn 


5. If the structural and node points are not identical, there are two options available 
in NPLDS for interpolating the aerodynamic forces from the aerodynamic nodes to 
the structural nodes. 

Option 1: (OPT1) 


Required data is as follows: 

1. Card input of the force weighting matrix [P] (weighting matrix of aerodynamic 
forces from aerodynamic nodes to structural nodes). 

2. FpL(q), FpL(q), Fpjjdg), read from the EOM tape 
The program calculates: 


[M 4 ] = SCALR2 [P] [Fp L (q>] 
[M 5 ] = SCALR2 [P] jF pL (q*)J 
[?] = SCALR3 [P] [F pL (a g )j 



r 


where the size of: 

[p] is the number structural loads by the number of aerodynamic panels 

[F PL (q)] is the number aerodynamic panels by the number of modes 

[F PL (q)] is the number aerodynamic panels by the number of modes 

[FpL(ag)] is the number aerodynamic panels by the number of gust penetration 

panels 

Limitations: 

Number of structural loads ^ structural nodes < 100 

Number of aero panels 100 number structural node 

Number of modes ^ 70 

Number of gust penetration panels ^ 35 

Number of structural nodes *£ 100 

Option 2: (OPT2) 

Required data is as follows: 

1. Read card input structural areas corresponding to structural nodes: [a s ]. 

2. Read structural (X, Y, Z)g coordinates from INTERP. 

3. Read local aero (X,Y)l coordinates from EOM. 

r " 

Ai 

4. Read local aero areas [A]l =' ' 'from EOM. 

V 

5. Read Fp^Cq), FpjJq), FpjJag) from EOM. 

Program operations are as follows: 

1. Calculate: 


[p] L = 
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where: 


Fin = FpL(q), F PL (q), F PL( a g ) 

[P]l= matrices of aerodynamic pressures at aerodynamic nodes for q, q, and dg. 

2. Call PLATEI using [P] L and (X,Y) L . 

3. Get SA from INTERP, extract transformation. 

4. Insert transformation in SA array from step 2. 

5. Call AINTG using (X,Y,Z) to interpolate [p] s (aerodynamic pressures on 
structural areas). 

6. Calculate aerodynamic force at structural nodes due to q, q, and a g as 

[FpL(q>] . [F PL (q)] , lFpL(«g)l = 

s s s 


\ 


\ 


\ 


V 


a j 


[p], 


where [P]s i s the aerodynamic pressure due to q, q\ and a g 
7. Then: 


[M 4 ] = SCALR2 [F PL (q)] s 


Similarly 


[M 5 ] = SCALR2 [ F pL (q)] 
[0] = SCALR3 [F pL (c? g )J s 


dp dp 

*Warning: This option is only valid over regions where and jj 1 - are equal or 

approximately equal to constants. y 



NPTAP and/or PTAP Final Output/LOAD-SET 


[Header] 


[M 3 1 

[m 4 ]1 

[M 5 ]J-k=l 

(h J 

[m 4 ]1 

[M 5 ]lk=2 

[?] J 


l Load 
f set 1 



EOF 

[Header] 

[M 3 ] 

lM4.ll 
[M 5 ]| k= 


l Load 
I set 2 


M- 


Surface 1 


where: [M3] = 


M/ 


Surface IS 
Surface 1 


[M 4 ] = 


Surface IS 


Similarily for [M3] and [0] 
and for all values of k 


Maximum size of 

[M 3 ] = [M 4 ] = [M S ] is 100 x 70 

[0] is 100 x (2) (35) 


A LOAD-SET is a result of processing all surfaces requested following the $NPLDS 
card, but before the next $ card. 

Note: The final output of PLDS is identical to that for NPLDS described above except 
that M3 matrix is omitted. 
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3.5 PRIMARY OVERLAY (L218,4,0) - VBMT 
Purpose of VBMT 

The L218 (LOADS) fourth primary overlay is named VBMT. VBMT reads the specific 
card input data to calculate shears, bending moments, and torsions and writes the load 
matrices on the magnetic file LTAP. 

Analytical Steps of VBMT 

VBMT performs its task in the following steps: 

1. VBMT initialize FETS for disk storage. 

2. The subroutine VBMTA is called to check all data for this load-set. 

3. The subroutine OPENMS is called to initialize a random access file. 

4. The subroutine VBMTA is called to read data for one load for this surface. 

5. The subroutine VBMTC is called to calculate M3 for this load. 

6. The subroutine VBMTD is called to calculate M4, M5, and for this load. 

7. If more loads or surfaces are to be processed, jump to step 4. 

8. The subroutine VBMTF is called to merge the matrices and write the final tape 
LTAP. 

9. If another load-set is to be read, jump to step 2. 

10. Close all FETS and return scratch files. 

11. Return control to (L218,0,0). 

The macro flow chart of this overlay is shown in figure 8. The subroutines called are 
displayed in table 5. 

I/O Devices of VBMT 

VBMT reads card input. Geometry and modal data is read from tape SATAP as provided 
by L215 (INTERP) (ref. 3). Equations of Motion data is read from tape EOMLOD as 
provided by L217 (EOM) (ref. 6). J-matrix data is read from tape JTAPE as provided by 
(L218,l,0), RGEN. 

Note: J-matrix data on JTAPE can also be used by NPLDS, since the mass matrix is 
the first record of each file. 

Regular and general print options control the printed output. Loads matrices are 
written on tape LTAP in a format acceptable to L219 (EQMOD) (ref. 4) and L221 
(TEV156) (ref. 5). 
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Figure 8. — Macro Flow Chart of Overlay ( L218,4,0 ) VBMT 
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Table 5. — Routines Called by VBMT 


OVERLAY (L2 18, 4,0) 
PROGRAM VBMT 


FETADD+ 
FETDEL+ 
OPENMS* 
RETURN? + 
VBMT A / 


VBMTA1 

VBMTA2 


VBMTA3 


VEMTC 


VBMTC1 

VBMTC2 


READTP+ 
/ VBMTC3 
READMS* 
WRITMS* 


LTOGT+ 

REACTP+ 

VBMTCU 


+ Indicates a routine in the DYLOFLEX alternate subroutine library. 
* Indicates a routine in the FORTRAN subroutine library. 
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Table 5. — ( Concluded ) 


VBMTD 


VEMTF 


VBMTD 1 
VBMTD2 
VBMTD 3 


READTP+ 

| VBMTD 1 
VBMTD4 


(_ WRITM S* 

’ FETDEL* 

V3MTF1 READMS* 

PETURNF+ 

WRTETP+ 


FEAETP+ 
J READMS* 
WPITMS* 
VBMTC3 


VBMTC4 
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3.5.1 VBMT PROGRAMMING SPECIFICATIONS 


To calculate [M 3 ] for a given surface number (IS), the inertia forces at the structural 
nodes and in local axis are: 

rm ' J Wx] ' ]' 

[ Fy J|ocal axis = Wy ' + tmftJ K]‘ fmiXil K] 

f F z:l = r™ii + KyjJ Re x ] + rm i x iJ pe v ] 

Mocal axis L x J L YJ 

where [F x> ], [Fy.], and [F z .] are the matrices of inertia forces in the x, y, and z 

directions. </> x , 4>y, <£z> ><£0 and <I>q are obtained from the INTERP tape (SATAP) 

u x u y z 

and mj, mjXj, m^, and m\z\ are obtained from the J-matrix (general card input or 
tape). These local axis forces are transformed to the inertia axis by: 


where: 


>x 


~ F x 

F y 

= SCALS [T] 

F y 

F z 

inertia 

F z 


axis 


local 

axis 


F 


x 


11 


In 


[F x ] 


n = number of modes 
i = number of nodes 


F 


x 


il 


in 


[F y ] and [F z ] are similar to [F x ]. [T] = the coordinate transformation matrix using 
0 X , 0y> 0 Z from the geometry data provided by INTERP. 

The matrix is formed using the Euler transformation triad [X] [Y] [z] or some other 
combination, where [X], [Y], [Z] are the individual axis transformation matrices 
making up the triad (see sec. 4.0, vol. I). 

SCALS = 1.0 or SCALE from card 6.9.3 (sec. 6.3, vol. I) for the specific structural 
panels. 
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I 


The forces in the inertia axis are summed at all required nodes to calculate shears at 
point A (load station or dummy load station). 


' = SCALE 1 Z; 

x A 1 

F Xj + 2 j 

SCALED; F v 

J x D ^ 

’ = SCALE 1 Z; 

y A 1 


SCALED; F y 

J x Dj 

’ = SCALE 1 Z; 

4 A 

F Xj +2 j 

SCALED; F v 

J X D: 


where: 

i = required structural node numbers, 
j = required dummy node numbers 
SCALE 1 is from card set 6.0 (sec. 6.3, vol. I) 
and 


X Dj 


= SCALE 1 


Fx i 

dummy node forces 

y °j 

Q 

> 

ll 

= SCALE) 


F Vi 

at dummy node j 

Z °j 

ll 

< 

s 

JP 

= SCALE 1 

•Zi 

F *i 



Note: IfSCALEDj= +2.0 (see vol. 1, sec. 6.3, card 6.9), 


F = M v = M v = A., is set = 0 

y d- x d- /j D- y 

j u j j 


If SCALED] = -2.0, 


F v = F_ = M„ is set = 0 
x Dj z Dj yDj 


and SCALED] = | -2.0 
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The moments at the structural nodes are obtained by [j] [<£] (local axis) as follows: 
[M x .] = rmjZjJ [<p y ] + E'mjyiJ [<p z ] + MxxjJ [ 4>q J + [-mjXjyj.] [<p e J - hnjXjZjJ [< $> e ] 

[M ] =-['nijZjJ [<p x ] + ['nijXjj [<p z ] + rmjXjyjJ [<t> g ] + MyyjJ ( <j> 9 ] + fmjyjZjJ l<p e ] 
y i x y z 

[M z .] ='frn i y i 4 [0 X ] - fnijXjJ [<t> y ] - fmjXjZjJ [<p d ] + fmjyjZjJ [<j> e ] + flzzjJ [</>q ] 
i x y z 


These moments are transformed to the inertia axis system by: 


My 

x i 

M y. 


M x 

-SCALS [T] 

My 

M, 

inertia 

m 7 

_ il 

axis 

L Z J 


local 

axis 


The bending moments at point A (load station or dummy load station) are now obtained 
by summing over all required nodes. 


M y = SCALE1 2 (A F + A 7 F v . + M y . ) + 2 MX: 
X A i y i i i y i i J J 

M v = SCALE 1 2 (A x F 7 . - A 7 F x . + M v .) + 2 MY: 
y a i i i i i y i j J 

M 7 = SCALE 1 2 (-A F y . - A y F v . + M 7 . ) + 2 MZ: 

Z A i y i x i x i y i z i j 1 


where: 


\ = BSi - BS a 


Ay, = BBLj- BBL a 


A =WLj-WL A 


MS j -SCALED j (A yDj F2 Dj BA ZDjFyDj tM XD _ ) 
MY j = SCALED j (A XDj F Z[)j -A ZDj F XDj + M yDj ) 


MZj -SCALEDj (-A F - A F + M ) 
J. J J J J 


A X D . - BS Dj - BS A 
A y D . =BBL d . -BBL a 

%. =WL Dj - WL A 

Dj is the j^ 1 dummy node 
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Note: See the previous note on SCALEDj. 

The dummy node forces and moments at each dummy load station are saved for future 
use. 

For load stations, the shears and moments are now transformed to the requested 
orientation. The angles 0 X , 0 y , and 0 Z used to calculate [T 1 ] for the transformation 
are obtained from card input (see card 6.7, sec. 6.3, vol. I). 


V x" 


rv 1 

v x 

V y 

= [T->] 

v y 

V z 

Point A 

v z 


requested axis ! 



and: 


M x 


M x 

M y 

= [T->] 

My 

_M Z _ 

A 

M Zj 


Point A 
inertia axis 


A, inertia axis 


Thus, we have M3 at point A: 


[M 3 ] 


V x 


v 

V 


Vx l 

x n 

v y 




V z 


• 

• 

M x 

= 

• 

* 

M„ 


• 


V 




M z 


M z, 

M z 

L r\ 

— 
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For a surface with two load points (A and B), and with all components (V x , V y , V z , M x , 
My, and M z ), the structure of [M3] as placed on LTAP would be: 


[M 3 ] 


V v 


X A 


V 


X B 


v 


y A 


v 


VB 


v 7 


Z A 


V 7 


Z B 


M v 


X A 

► 

M v 


X B 


M , 


^A 


M , 


VB 


m 7 


Z A 


M, 


Z B 

— * " 

► 

— _ 



Surface I 


Surface N 


If one or more component is missing, the matrix would close up; also, each surface is 
treated separately and is merged onto the preceding surface. 


To calculate [M4], [M5], and {</>}, obtain the aero force coefficient matrices Fp^Cq), 
FpL(q), FpjJdg) from the Equations of Motion tape and transform the forces into the 
matrix axis system. Then: 


F x (q,q,a g ) 


r 1 

F y (q.q.<%) 

= SCALAj « 

Tj 

F z . (q,q,d„) 

inertia 

k. J 

1 * 

axis 



where {Tj} is the transpose of the third row of [T] with 0 Z 


local 

axis 


= 0, or: 


F x . (q,q,<ig) = SCALAj T(3, 1 ) F pL . (q,q,a g ) 
F y . (q,q,a g ) = - SCALAj T(3,2) F pL . (q,q,dfg) 
F z . (q,q,d g ) = SCALAj T(3,3) F pL . (q,q ,d g ) 


SCALAj - 1 .0 or SCALE from card 6.93 (section 6.3, volume 1 ) 
for the specific aerodynamic panels. 
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The aero forces in the inertia axis are summed at all required nodes to calculate shears 
at point A (load station or dummy load station, inertia axis): 


V Y = SCALE SF y + 2 SCALED; F Y 

X A i i j J x Dj 

V v = SCALE 2 F v + 2 SCALED; F v 

^A i Yi j • Vd j 

V 7 = SCALE 2 F 7 + 2 SCALED: F 7 

Z A i z i j J z Dj 


where: 

i = required aerodynamic node numbers 
j = required dummy node numbers 
SCALE = SCALE2 for M 4 and M 5 
= SCALE3 for Ifr 

dummy node forces 
at dummy node j 


x °j 

11 

< 

X 

J? 

= SCALEZj 

F x 

y Dj 


= SCALED; 

F v 

Z D: 

II 

< 

N 

U 

= SCALER 

Fz i 


Note: See the previous note on SCALED} . 


For slender bodies (Z bodies) where the force is in the Z direction: 


0 


inertia axis 


^inertia axis 

F = F 

z inertia axis z aero local axis 


For slender bodies (Y bodies) where the 

F v 


inertia axis 

^inertia axis 
F„ 


inertia axis 


force is in the Y direction: 

= 0 

y aero local axis 

= 0 
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The force coefficient matrices from EOM are partitioned for slender bodies to contain F y 
and F z . For example, Fp^Cq) for a slender body with y and z forces would appear as 
follows: 


[f P l (q)] = 


Fp L (Q)y 


F PL (q)z 


Note: Fy or F z may be zero, but both F y and F z will be on the magnetic File. 

This is possible since the maximum possible number of nodes for slender bodies is less 
than 50. Thus, the maximum number of nodes for [F y and F z ] would be less than 100. 

The bending moments at point A (load station or dummy load station) are obtained by 
summing over all required nodes: 

M = SCALE Z (A v . F z . + A 7 F v ) + I MX: 

A i >i z i z i V] j J 

. M yA = SCALE Z (A x . F 2 . - A z . F x .) + Z M Yj 

M = SCALE 2 (-A F x . - A X .)F V + 2 MZ: 

a i yj ; j 


where: 


^ Xi =BS i- BS A 

Ay. = BBLj - BBL^ 

A z . = WLj - WL A 

M Xj ” SCALEDj (AyDj Fzp>. + ^zjj. ^yj). + ^xj). ^ 
M yj = SCALEDj (A X[) ^ F zd . - A z ^ F x ^ + M yo . ) 
M Z j = SCALEDj (-Ay D _ F x ^ - A XJ) _ Fy^ + M z ^) 
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A X D . BS Dj - BS A 




— BBLq. - BBL^ 


J 


\ ' WL Dj 


- WL, 


Dj is the dummy node 

SCALE = SCALE2 for M 4 and M 5 
= SCALE3 for 0 

Note: See the previous note on SCALEDj. 

The dummy node forces and moments at each dummy load station are saved for future 
use. For load stations, the shears and moments are then transformed to the requested 
orientation (for a dummy load there is no transformation). The angles 0 X , 6 y , and 6 Z are 
obtained from card input (see card set 6.0, sec. 6.3, vol. I). 


and: 


V x" 


V x 

v y 

= [T-‘ ] 

v y 

v z 

A, requested 

_ V Z _ 

L_ — 

axis 


Mx 


M X “ 

M y 

= [IT 1 ] 

My 

_M Z _ 

A, requested 
axis 

M x 


A, inertia 
axis 


A, inertia 
axis 


Thus: 



f V xl 

V y 


pv v 1 

x x n 



= 

• 

II 

< 

S 

M x 

M y 


• 


M z 

A 

m 7 m 7 

Z 1 z n 


where n = number of modes. Similarly for M 5 : 


[0] 


where g = number of gust zones. 


V v 


V 

V 

x 

v y 


X 1 

x g 

V z 

= 

# 


M x 

M_ 


M 

M 

z 

A 

M Zj 

M 

gj 
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VBMT Output TAPE (LTAP) 

Note: A load set is the result of processing all surfaces requested following the 
LOAD-SET card, but before the next LOAD- SET or $ card. 

[ Header] 

[M 3 ] 

[M 4 ] 


[M 5 ] 

[£ ] 


} 


k=l 


[M 4 ] 

[M s ] f k=Nk 
[% 1 
EOF 
[Header] 

[M3] 
t“4l , 


r Load set 1 


• j 


r Load set 2 


Maximum size of [M 3 ] = [M 4 ] = [M 5 ] is 100 x 70; of [<£] is 100 x 2(35). 
Expansion of [M 3 ] on LTAP (the first output on LTAP): 


[M 3 ] = 


Vx l 

point 1 

V x 

point m 

m 


• 

. ► 


point 1 

m 7 

point m 


point 1 

. 

• ► 

M Zn 

point n 


; 1 

• 

1 


Y Surface 2 
Surface 1 

Similarly for [M 4 ], [M 5 ], and [<£] and for all values of k. 


r Surface 1 
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Local to Reference Coordinates 


The coordinate locations (X,Y,Z) of the loads from card 6.9.2 (sec. 6.3, vol. I) are 
always in reference coordinates. Thus it can be directly compared with the 
reference coordinates from INTERP to determine the structural nodes to be 
included in the summation for [M3]. 

However, the coordinates (X,Y) associated with the aero data (from the EOM 
tape) are in local coordinates. They must be changed back to structural 
coordinates before being compared to load coordinate _locations _QC,Y,Z) in the 
selection of the aero nodes to get summation for [M4], [M5], and [(f)]. 

This is accomplished by LTOGT from the DYLOFLEX library. CALL LTOGT 
(X,Y,Z, NPTS, R,T), where (X,Y,Z) input is (X,Y) from EOM with Z being zero, 
and R and T are from the SA array from INTERP (rotation and transformation). 
LTOG will then return (X,Y,Z) in structural coordinates. (LTOGT is an entry 
point in AINTT). 

For slender bodies, the Z coordinate in (X,Y,Z) from LTOGT may be replaced by 
the Z from input card 6.6. 


3.6 DATA BASES 

L218 (LOADS) data bases include input and output files plus internal scratch 
files and common block storage. 

3.6.1 INPUT DATA 

The input data is from two sources, cards and magnetic files. 

Card Input Data 

For a complete description of all the card input formats, see section 6.3 in 
volume I of this document. 

Tape Input Data 

For a complete description of the tape input data see section 6.4 in volume I of 
this document. 

3.6.2 OUTPUT DATA 

The output data may be of two types, printed and magnetic files. 

Printed Output Data 

For a complete description of the printed output data, see section 6.5.1 in 
volume I of this document. 
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Magnetic Files (Tape or Disk) 

For a complete description of the magnetic file output data, see section 6.5.2 in 
volume I of this document. 

3.6.3 INTERNAL DATA 

Common blocks and blank common (dynamic storage) are used to pass data from 
one routine to another within an overlay. Temporary (scratch) disk files are used 
in (L218,2,0), AVD. Random access disk files are used in (L218,3,0), NPLDS, and 
(L218,4,0), VBMT. 

AVD Internal (Temporary) Disk Storage 

FETADD is used to initialize buffer storage. FETDEL and RETURNF are called 
to delete the files after they are no longer needed. Three files are used: 

o IFMj contains merged data. 

o IFM2 contains merged M2 data. 

o IFM3 contains merged M3 data. 

These data are written on these three files in subroutine MERGE with 
FORTRAN write statements and read in subroutine AVDTAP. The files are 
initialized and deleted in AVD. The record structure (fig. 9) for IFMj, IFM2, and 
IFM3 is identical. 

NPLDS/PLDS Internal (Temporary) Disk Storage 

All temporary disk storage in this overlay is accomplished on the file MERGMB 
using random access methods. The file MERGMB is initialized and deleted in 
NPLDS. The subroutines that call WRITEMS to write on MERGMB are NPLDB2 
and NPLDD4- NPLDB2 writes the matrix M3, and NPLDD4 writes the matrice 
M4, M5, and (j>. The subroutine that calls READMS from MERGMB is NPLDH1. 

The matrices or records are indexed with the following numbers: 

Matrix Random access key 


Mg 

1 

to ISMAX 

M4 

ISMAX + 1 

to 21(ISMAX) 

m 5 

2KISMAX) + 

1 to 41(ISMAX) 

4> 

4HISMAX) + 

1 to 61(ISMAX) 


where ISMAX is the maximum number of surfaces. 

The read/write activity on MERGMB in NPLDS/PLDS is displayed in figure 10. 


Size 


JftU 


NRCWS, 
NROWS f 


X MODES 
X MODES 


NRCWS^ X MODES 


NROWS X MODES 
NSURF 



where [M 3 ] is the [M 3 ] for surface i 
NSURF is the last surface 
MODES is the number of modes 

NROWS is dependent on the number of selected nodes 
and the selected matrix components. 
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Figure 9. — Record Structure of fFM p IFM 2 , and IFMj 



Record No. 

Record 

Size 

Written In 

Read In 

Contents 

1 

NODES ( 

MXMODE+1) 

NPLE 

)B2 

NPL 

DH1 

Mj - Surface 1 

ISMAX 







iT 3 - Surface ISMAX 

ISMAX+1 



NPLE 

)D4 



ff 4 - Surface 1, freq 1 

+2 







freq 2 
• • 

: +20 







• freq 20 

+21 







Surface 2, freq 1 
• 

21 ( ISMAX) +1 







• 

M 5 - Surface 1, freq 1 
• 

41 (ISMAX) +1 

f 

NODES (COL+1) 

i 





__ • 

$ - Surface 1, freq 1 

• • 

61 (ISMAX) 

i 






1" - Surface ISMAX, freq 20 


Figure 10. — Read /Write Activity on MERGMB in NPLDS/PLDS 











VBMT Internal (Temporary) Disk Storage 

All temporary disk storage in this overlay is accomplished on the file MERGMB 
using random access methods. MERGMB is initialized and deleted in VBMT. The 
subroutines that call WRITEMS to write on MERGMB are VBMTC2, VBMTD3, 
and VBMTD4. The subroutines that call READMS to read from MERGMB are 
VBMTC2, VBMTD4, and VBMTF1. The read/write activity on the file MERGMB 
in VBMT is displayed in figure 11. 

Common Blocks 

Table 6 displays the common blocks used in the program and the overlays in 
which they are used. 

The LABELED common blocks are used for communication between the main and 
primary overlays, and for communication between routines in a primary overlay. 
The block names and contents are described in table 7. The "T” heading in table 7 
refers to variable type. The codes used are as follows: 

I Integer 

R Real 

C Complex 

L Logical 

H Hollerith 

Blank common is used in the primary overlays as a variable length working 
storage area. The length of required arrays is calculated, and the first word 
address and variable dimension of the array is passed through the subroutine 
calling sequence for those routines which require it. 
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Record No. 

Record Size 

Written In 

Read In 

Contents 

1 

MXLOAD 
MXLOAD+ 1 

MXLOAD (1+NK) +1 

MXLOAD (2+NK)+l 

MXLOAD (3+NK) +1 

• 

• 

LAST 

I ROW (MX! 

MODE+2) 

i 

VBM 

TC2 

NPLI 

1 

)Hl | 

M 3 for load No . 1 

M h for load No. 1, freq 1 

M 5 for load No. 1, freq 1 
• 

• 

-A. 

$ for load No. 1, freq 1 

• 

• 

Dummy load No . 1 

• 

• 

• 


where 

MXLOAD 

NK 

I ROW 

MXMODE 

LAST 

MXDUM 


is the number of loads for this load set 

is the number of frequencies 

is the number of rows in this matrix 

is the number of modes 

= MXLOAD ( 3+NK) +MXDUM ( 1+NK) 

is maximum number of dummy loads 


Figure 11. — Read /Write Activity on MERGMB in VBMT 




Table 6. — Common Blocks in Each Overlay 


V XC0MM0N 

X.BLOCK 

OVERLAY'S. 

CGEN 

CGEN1 

CGFN2 

RWBUFF 

CAVD 

CLOC 
CITMRM 
CMAT 1 

CARRAY 
Q3 

CNPLP1 
CNPLD2 
CVBMT1 
C VBMT 2 
CVBMTF1 

12 

1-3 

CQ 

(L218 ,0,0) 
L218 

XX X 



M 

(L218, 1,0) 
RGEN 

X X X X 



\r 

dS. 

(L218, 2,0) 
AVD 

X X X X X 

X X X X X 


X 

(L218, 3,0) 
NPLDS 

X X X X 


X X 

X 

(L218 ,4,0) 
VBMT 

XX X 


XXX 

\r 
























Table 7. — Contents of Common Blocks 


BLANK COMMON: 

AVD Dynamic Storage 

DESCRIPTION : 

Variably 

dimensioned arrays 









_ _J 

NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

INODE 

I 

100 

i 

Contain selected nodes 

2 

LFMl 

I 

LBUF 


FET buffer for IFM1 

3 

LFM2 

I 

LBUF 


FET buffer for IFM2 

4 

LFM3 

I 

LBUF 


FET buffer for IFM3 

5 

LARRAY 

I 

LBUF 


FET buffer for IDISK 

6 

XYZCOR 

R 

3*NODES 

BS ,BBL V , 
^X 

( X , Y , Z ) coordinates for sur- 
face IS (card input) 

7 

ALB 

R 

NODES 

LB 

(LB) interpolation coefficient 

8 

ALT 

R 

NODES 

LT 

(LT) interpolation coefficient 

9 

ALTT 

R 

NODES 

LTT 

(LTT) interpolation coefficient 

10 

ISUB 

I 

NODES 


Requested nodes 

11 

GEOM 

R 

3 *NODES 

BS ,BBL t , 

WLj 

( X , Y , Z ) coordinates from 
INTERP 

12 

BS2 

R 

NODES 

BS i + l 

X coordinate from INTERP for 
the (1+1) node 

13 

THXYZ 

R 

3*NODES 

9 x' e y' e Z 

0 x ,6 y' e z ^° r sur ^ ace IS 

14 

PX 

R 

NODES* 

IMXMODE 

4 

x i 

<P 

Y x 

15 

PX2 

R 

NODES* 

MXMODE 

<f> 

x i+l 

<f> x for (1+1) node 

16 

PY 

R 

NODES* 

MXMODE 

<J> 

*i 

<t> 

y 

17 

PY2 

R 

NODES* 

MXMODE 

y i + l 

for (1+1) node 

18 

PZ 

R 

NODES* 

MXMODE 

♦z. 

1 

*z 

19 

PZ 2 

R 

NODES* 

MXMODE 

+ z 

Z i+1 

^z for (1+1) node 

20 

PPX 

R 

NODES* 

MXMODE 

CD 

H* 

^ex 
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VARIABLE 
21 PPX2 
2 2 PPY 

23 PPY2 

24 PPZ 

25 PPZ 2 

26 SA 


DIM. 


NODES* 

MXMODE 


NODES* 

MXMODE 

*e yi 

NODES* 

MXMODE 

* e A+i 

iNODES* 

MXMODE 

* ez i 

NODES* 

!mxmode 

* 6z i+i 


DESCRIPTION 


9 6x for (1+1) node 
^0y 

^0y for (1 + 1) node 
^ 0 z 

^ 0 z for (1 + 1) node 
SA array from INTERP 







Table 7. — (Continued) 


LABELED COMMON NAME: CGEN 

DESCRIPTION: 

Contains 

general 

data required bv all overlays. 















NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

TDUM 

H 

8 



Temporary input array 

2 

INFIL 

I 

1 



Name of card input file (=5) 

3 

IUTFIL 

I 

1 



Name of card output file (=6) 

4 

IGPRNT 

I 

1 



General print control 

5 

KMOD 

I 

1 



Module option control (AVD=1, 
NPLDS=2 , PLDS=3, VBMT=4) 

6 

I OUT 

I 

1 



Name of temporary general output 
file . 

7 

IGS 

I 

1 



Load-Set counter 

8 

NERROR 

I 

1 



Fatal Error Counter 

9 

ISMAX 

I 

1 



Maximum Number of surfaces for 
this load-set 

10 

WAR 

H 

3 



Warning Message 

11 

FAT 

H 

3 



Fatal Error Message 

12 

LBUF 

I 

i 



Dimension value for OUTBUF-buf- 
fer length. 

13 

IAVDl 

H 

i 



AVD output tape name 

14 

NPTAP 

H 

i 



tfPLDS output tape name 

15 

IPTAP 

H 

i 



PLDS output tape name 

16 

LTAP 

H 

i 



^BMT output tape name 

17 

IEOMLD 

H 

i 



File name for EOM tape from 
L217 (EOM) 

18 

IDISK 

H 

i 



File name for interpolation 
tape from L215(INTERP) 

19 

MASS TP 

H 

i 



File name for mass tape 

20 

JTAPE 

H 

i 



File name for the J-Matrix 
tape 

21 

IXYZ 

I 

i 



Lndicator defining transforma- 
tion order 

22 

LC 

I 

i 



Line Count for FETADD 
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LABELED COMMON NAME: CAVD 


DESCRIPTION: AVD subroutine options and scale factors 


NO J 


VARIABLE 


DIM. 


ENG 


I CARD 


IIS 


3 

4 

5 

6 

7 

8 
9 


12 

L3 

L4 


INTER 
IPRINT 
I ROT 
IS 

IUNIT 

NIS 

NRR 

SCALR1 

SCALR2 

SCALR3 

IANGLE 

ISP 


I 

I 

I 

I 

I 

I 

I 

R 

R 

R 

I 

I 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


. NOM 


DESCRIPTION 


=1 for matrix input by card; 

=0 otherwise 

=1 if another surface is to be 
input; =0 otherwise 

Interpolation control 

Print control 

=1 for rotation; =0 otherwise 
Current surface number 
=1 for Metric; =2 for English 
Not used 

Diagnostic error variable 
Matrix scale factor 
Matrix scale factor 
Matrix scale factor 
Theta input control 
Previous surface number 









LABELED COMMON NAME: 


CGEN1 















Table 7. — (Continued) 


LABELED COMMON NAME: CITMRN 


DESCRIPTION: AVD subroutine translation and rotation direc- 

tives . 


NO. 


VARIABLE 


DIM. 


ENG. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 
17 


ITM3X 

ITM2X 

ITMlX 

ITM3Y 

ITM2Y 

ITM1Y 

ITM3Z 

ITM2Z 

ITMlZ 

IRM3X 

IRM2X 

IRM1X 

IRM3Y 

IRM2Y 

IRM1Y 

IRM3Z 

IRM2Z 

IRM1Z 



NOMJ 


DESCRIPTION 


Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 

Code 

for 


translational matrix 
translational matrix 
translational matrix 
translational matrix 
translational matrix 
translational matrix 
translational matrix 
translational matrix 
translational matrix 
rotational matrix 
rotational matrix 
rotational matrix 
rotational matrix 
rotational matrix 
rotational matrix 
rotational matrix 
rotational matrix 
rotational matrix 



LABELED COMMON NAME: CLOC 

DESCRIPTION : AVD subroutine matrix sizes and locations 


NO. 

VARIABLE 

D 

DIM. 

ENG. NOM. 

i 

NODES 

I 

i 


2 

MXNODE 

I 

i 


3 

MXMODE 

I 

i 


4 

LXYZ 

I 

1 


5 

LLB 

I 

i 


6 

LLT 

I 

i 


7 

LLTT 

I 

i 


8 

LCODE 

I 

1 


9 

LISUB 

I 

1 


10 

LGEOM 

I 

1 


11 

LBS 2 

I 

1 


12 

LTH 

I 

1 


13 

LPX 

I 

1 


14 

LPX2 

I 

1 


15 

LPY 

I 

1 


16 

LPY2 

I 

1 


17 

LPZ 

I 

1 



DESCRIPTION 


Number of nodes 
Total number of 
Total number of 

requested 

nodes 

modes 

Dynamic 

XYZCOR 

storage 

location 

for 

Dynamic 

ALB 

storage 

location 

for 

Dynamic 

ALT 

storage 

location 

for 

Dynamic storage 
ALTT 

Not Used 

location 

for 

Dynamic 
I SUB 

storage 

location 

for 

Dynamic 

GEOM 

storage 

location 

for 

Dynamic 

BS2 

storage 

location 

for 

Dynamic 

THXYZ 

storage 

location 

for 

Dynamic 

PX 

storage 

location 

for 

Dynamic 

PX2 

storage 

location 

for 

Dynamic 

PY 

storage 

location 

for 

Dynamic 

PY2 

storage 

location 

for 

Dynamic 

PZ 

storage 

location 

for 


13 


LPZ2 


I 


1 


Dynamic storage location for 
PZ2 










Table 7. — ( Continued ) 


LABELED COMMON NAME: CLOC (continued) 

DESCRIPTION : ___ 


NO. 

VARIABLE 

19 

LLPX 

20 

LLPX2 

21 

LLPY 

22 

LLPY2 

23 

LLP 2 

24 

LLPZ2 

25 

LTM3X 

26 

LTM2X 

27 

LTM1X 

28 

LTM3Y 

29 

LTM2Y 

30 

LTM1Y 

31 

LTM3Z 

32 

LTM2Z 

33 

LTM1Z 

34 

LRM3X 

35 

LRM2X 

36 

LRM1X 

37 

LRM3Y 

38 

LRM2Y 

39 

LRM1Y 

40 

LRM3Z 

41 

LRM2Z 

42 

LRM1Z 




DESCRIPTION 




Dynamic storage 
PPX 

Dynamic storage 
PPX2 

Dynamic storage 
PPY 

Dynamic storage 
PPY 2 

Dynamic storage 
PPZ 

Dynamic storage 
PPZ 2 

Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not used 
Not ‘used 
Not used 
Not used 
Not used 
Not used 


location for 
location for 
location for 
location for 
location for 
location for 



















Table 7. — (Continued) 


(labeled COMMON NAME: CMAT 



DESCRIPTION: 

AVD subroutine 

problem size. 















2 
1 O 

9 

VARIABLE 

T 

DIM. 

ENG. 

NOM. 

DESCRIPTION 

1 

2 

3 

4 

5 

6 

ITOTAL 

ITOTIS 

ISTOT 

ISM3 

ISM2 

ISM1 

I 

I 

I 

I 

I 

I 

i 

1 

1 

1 

100 

100 

100 



Total number of rows for this 
load-set . 

Total number of rows for sur- 
face IS. 

Current number of surfaces for 
this load-set. 

Total rows written on FM3 
per surface 

Total rows written on FM2 
per surface 

Total rows written on FM1 
per surface 


15 



Table 7. — (Continued) 


LABELED COMMON NAME: CNPLDl 


DESCRIPTION : 

NPLDS/PLDS subroutine array sizes, location and 



scale factors. 








NO. 

VARIABLE 

T 

DIM. 

ENG. NOMj 

DESCRIPTION 

1 

NODES 

I 

1 


Number of nodes requested 

2 

KXNODE 

I 

i 


Total number of nodes 

3 

MXMODE 

I 

1 


Total number of modes 

4 

LNM3 

I 

1 


Dynamic storage location for 
NM3 

5 

LNM4 

I 

i 


Dynamic storage location for 
NM4 

6 

LMM5 

I 

1 


Dynamic storage location for 
MM 5 

7 

LNC3 

I 

1 


Dynamic storage location for 
NC3 

3 

LIEOM 

I 

i 


Buffer location for IEOMLD 

9 

LIDISK 

I 

1 


Buffer location for IDISK 

10 

LNPTAP 

I 

1 


Buffer location for NPTAP 

11 

LMERG 

I 

1 


Buffer location for OPENMS 

12 

LMASS 

I 

i 


Buffer location for MASSTP 

13 

LPTAP 

I 

i 


Buffer location for IPTAP 

14 

LINODE 

I 

i 


Dynamic storage location for 
INODE 

15 

LXMASS 

1 I 

i 

1 

Dynamic storage location for 
XMASS 

16 

LAREA 

I 

i 


Dynamic storage location for 
AREA 

17 

LPMAT 

I 

i 


Dynamic storage location for 
PMAT 

18 

LGEOM 

I 

i 


Dynamic storage location for 
! GEOM 

19 

LPZ 

I 

i 


; Dynamic storage location for 
1 PZ 

20 

LXYL 

I 

i 


Dynamic storage location for 
XYL 

- - 
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Table 7 . — (Continued) 


LABELED COMMON NAME: 

CNPLDl 

(continued) 

DESCRIPTION : 

















NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

21 

LAREAL 

I 

1 


Dynamic storage location for 
AREAL 

22 

LXKVAL 

I 

1 


Dynamic storage location for 
XKVAL 

23 

LFPL 

I 

i 


Dynamic storage location for 
FPL 

24 

LSA 

I 

1 


Dynamic storage location for 
SA 

25 

LAST 

I 

1 


Last location of dynamic stor- 
age 

26 

NK 

I 

1 

NK 

Number of frequencies 

27 

NAERO 

I 

1 


Number of aeropanels 

28 

NLOAD 

I 

i 


Number of structural loads 

29 

ICARD 

I 

1 


=1 for matrix input by cards; 
-0 for matrix input by tape. 

30 

IIS 

l 

I 

1 


=1 for another surface; =0 
otherwise 

31 

I OPT 

I 

1 


Option control 

32 

IPRINT 

I 

1 


Print control 

33 

IS 

I 

1 


Current surface number 

34 

ISP 

I 

1 


Previous surface number 

35 

ISTOT 

I 

1 


Number of surfaces for this 
load-set 

36 

j 

ITOTAL 

I 

1 


Total number of rows for this 
load-set 

37 

IUNIT 

I 

1 


=1 for Metric; =2 for English 

38 

IM3 

I 

1 


Random Access Key Number - Ms 

39 

IM4 

I 

1 


Random Access Key Number - M 4 

40 

IM5 

I 

1 


Random Access Key Number - M 5 

41 

IC3 

I 

1 


Random Access Key Number - 

42 , 

NGUST 

I 

1 

g 

Number of gust panels 





NO. VARIABLE 


4 3 LROW 

44 MERGMB 

45 SCALR1 

46 SCALR2 

47 SCALR3 

48 UNIT 


Dynamic storage location for 
NROW 

Random Access File 
SCALR1 Scale factor for M 3 

SCALR2 Scale factor for M 4 and M s 

SCALR3 Scale factor for 

Metric or English units 
















Table 7. — (Continued) 


LABELED COMMON NAME: CVBMT1 

DESCRIPTION: 

VBMT subroutine option, array sizes, and 



locations . 







NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

i 

ICOMP 

I 

i 


Current component number 

2 

IC3 

I 

1 


Random Access Key - <T 

3 

IIS 

I 

1 


Type of input indicator 

4 

IL 

I 

1 


Current sequencial load number 

5 

ILD 

I 

i 


Current sequencial dummy load 
number 

6 

ILDS 

I 

i 


Number of Dummy Nodes to sum to 
tnis surface 

7 

ILDT 

I 

i 


= 0 for Dummy Load; = 1 otherwise 

8 

ILL 

I 

i 


Current load number (not se- 
guencial 

9 

ILS 

I 

1 


Load-Set number 

10 

IMD 

I 

i 


Random Access Key - Dummy Loads 

11 

IMXX 

I 

i 


=1 for MX comp; =0 otherwise 

12 

IMYY 

I 

i 


=1 for MY comp; =0 otherwise 

13 

IMZZ 

I 

i 


=1 for MZ comp; =0 otherwise 

14 

IM3 

I 

i 


Random Access Key - M 3 

15 

IM4 

I 

i 


Random Access Key - 

16 

INTMP 

I 

i 


Temporary input file 

17 

IMS 

I 

i 


Random Access Key - Ml 

18 

IPRINT 

I 

i 


Print control 

19 

IR 

I 

i 


Error indicator 

20 

IRR 

I 

i 

i 

Error indicator 

21 

IS 

; I 

i 


Current surface number 

22 

ISLEND 

I 

i 

1 

=1 for a slender body; =0 
otherwise 

23 

ISP 

I 

i 


Previous surface number 

24 

ISTOT 

I 

i 


Number surfaces for this load- 
set 
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Table 7. — (Continued) 

LABELED COMMON NAME: CVBMT1 (continued) 

DESCRIPTION: _ 


NO. 

VARIABLE 

25 

ITOTAL 

26 

IUNIT 

27 

IVBMTA 

28 

IVX 

29 

IVY 

30 

IVZ 

31 

IZ 

32 

LAST 

33 

! LAXYZ 

34 

LDAX 

35 

LDAY 

36 

: LDAZ 

37 

LFPL 

38 

LGEOM 

39 

LICOM 

40 

LIDISK 

41 

LIDNOD 

42 

LIEOM 

43 

LI I Z 




DESCRIPTION 


Total number rows on final 
matrix 

=1 for Metric; =2 for English 

=1 to check data for load-set 
=2 to read a load 

= 1 for VX computed; =0 otherwise 

= 1 for VY computed; =0 otherwis 

= 1 for VZ computed; =0 otherwis 

Number Z override coordinate 

Last location of dynamic stor- 
age 

Dynamic storage location for 
AXYZ 

Dynamic storage location for 
DAX 

Dynamic storage location for 
DAY 

Dynamic storage location for 
DAZ 

Dynamic storage location for 
FPL 

Dynamic storage location for 
GEOM 

Dynamic storage location for 
ICOM 

Buffer storage for IDISK 

Dynamic storage location for 
IDNOD 

Buffer storage for IEOMLD 

Dynamic storage location for 
IZ 



Table 7. — (Continued) 


LABELED COMMON NAME: CVBMT1 (continued) 

DESCRIPTION: 















NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 


44 

LILSEQ 

I 

i 


Dynamic storage location 
ILSEQ 

for 

45 

LINODA 

I 

1 


Dynamic storage location 
INODA 

for 

46 

LINODS 

I 

1 


Dynamic storage location 
INODS 

for 

47 

LINTMP 

I 

i 


Buffer storage location for 
INTMP 

48 

LISEQ 

I 

1 


Dynamic storage location 
ISEQ 

for 

49 

LISLOD 

I 

i 


Dynamic storage location 
IS LOAD 

for 

50 

LJTAPE 

I 

1 


Dynamic storage location 
JTAPE 3uf f er 

for 

51 

LLNODD 

I 

1 


Dynamic storage location 
LNODD 

for 

52 

LLNODE 

I 

i 


Dynamic storage location 
LNODE 

for 

53 

LLTAP 

i 

I 

1 


Dynamic storage location 
LTAP Buffer 

for 

54 

LMASST 

I 

1 


Dynamic storage location 
MASS TP Buffer 

for 

55 

LMERG 

I 

1 


Dynamic storage location 
MERGMB Index 

for 

56 

LMMXX 

I 

1 


Dynamic storage location 
MMXX 

for 

57 

LMMYY 

I 

1 


Dynamic storage location 
MMYY 

for 

58 

LMMZZ 

I 

1 


Dynamic storage location 
MMZZ 

for 

59 

LMVX 

I 

1 


Dynamic storage location 
MVX 

for 

60 

LMVY 

I 

1 


Dynamic storage location 
MVY 

for 
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Table 7. — (Continued) 


LABELED COMMON NAME: CVBMT1 (continued) 

DESCRIPTION: 

















NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

61 

LMVZ 

I 

i 


Dynamic storage location for 
MVZ 

62 

LUC 3 

I 

i 


Dynamic storage location for 
NC3 

63 

LNMD 

I 

i 


Dynamic storage location for 
NMD 

64 

LNM3 

I 

1 


Dynamic storage location for 
NM3 

65 

LNM4 

I 

i 


Dynamic storage location for 
NM4 

66 

LNM5 

I 

1 


Dynamic storage location for 
NM5 

67 

LNODEZ 

I 

1 


Dynamic storage location for 
NODEZ 

68 

LPPX 

I 

i 


Dynamic storage location for 
PPX 

69 

LPPY 

I 

J- 


Dynamic storage location for 
PPY 

70 

LPPZ 

I 

i 


Dynamic storage location for 
PPZ 

71 

! LPX 

I 

i 


Dynamic storage location for 
PX 

72 

LPY 

I 

i 


Dynamic storage location for 
PY 

73 

LPZ 

I 

i 


Dynamic storage location for 
PZ 

74 

LROW 

I 

! 

i 


Dynamic storage location for 
NROW 

75 

LSAX 

I 1 

i 

1 

Dynamic storage location for 
SAX 

76 

LSAY 

i 

i 


Dynamic storage location for 
SAY 

77 

LSAZ 

I 

i 


Dynamic storage location for 
SAZ 
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Table 7. — (Continued) 


LABELED COMMON NAME: CVBMT1 (continued) 

DESCRIPTION: 

















NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

78 

LSCALA 

I 

i 


Dynamic storage location for 
SCALA 

79 

LSCALD 

I 

1 


Dynamic storage location for 
SCALD 

80 

LSCALS 

I 

1 


Dynamic storage location for 
SCALS 

31 

LTHEX 

I 

1 


Dynamic storage location for 
THEX 

32 

LTHEY 

I 

i 


Dynamic storage location for 
THEY 

83 

LTHEZ 

I 

1 


Dynamic storage location for 
THEZ 

84 

LXIXX 

I 

1 


Dynamic storage location for 
XI XX 

85 

LXIYY 

I 

i 


Dynamic storage location for 
XIYY 

86 

LXIZZ 

I 

i 


Dynamic storage location for 
XIZZ 

87 

LXKVAL 

I 

i 


Dynamic storage location for 
XKVAL 

88 

LXMA 

I 

i 


Dynamic storage location for 
XMA 

89 

LXMX 

I 

i 


Dynamic storage location for 
XMX 

90 

LXMXY 

I 

i 


Dynamic storage location for 
XMXY 

91 

LXMXZ 

I 

i 


Dynamic storage location for 
XMXZ 

92 

LXMY 

I 

i 


Dynamic storage location for 
XMY 

93 

LXMYZ 

I 

i 


Dynamic storage location for 
XMYZ 

94 

LXMZ 

I 

i 


Dynamic storage location for 
XMZ 
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Table 7. — ( Continued ) 


LABELED COMMON NAME: CVBMT1 


DESCRIPTION: 


— 

NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

95 

LZ 

I 

1 


Dynamic storage location for 

n 

Li 

96 

MERGMB 

I 

1 


Random Access File 

97 

MX AERO 

I 

i 


Maximum nodes from EOM 

93 

MXCOMP 

I 

1 


Maximum Load Components for 
this surface 

99 

MXDUM 

I 

1 


Maximum number of dummy loads 

100 

MXLOAD 

I 

1 


Number of loads for this load- 
set 

101 

MXMODE 

I 

1 


Number of modes 

102 

MXNODE 

I 

1 


Maximum number of nodes from 
INTERP 

103 

NAERO 

I 

i 


Number of aero nodes requested 

104 

NGUST 

I 

i 


Number of gust panels 

105 

NK 

I 

1 


Number of frequencies 

106 

NODES 

1 

I 

1 


Number of structural nodes 
requested 

107 

SCALE1 

R 

1 

SCALE 1 

Scale factor for M^ 

108 

SCALE2 

R 

1 

SCALE 2 

Scale factor for FT 4 and M” 5 

109 

SCALE3 

R 

1 

SCALE 3 

Scale factor for $ 

110 

T 

R 

(3/4) 

i 

i 

i 

Transformation and rotation 
matrix from SA array 
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Table 7. —(Continued) 


LABELED COMMON NAME: CVBMTFl 


DESCRIPTION: 


Maximum, first and last surface numbers 


NO. VARIABLE 



ENG . NOM . 


DESCRIPTION 


Maximum surface number in ISEQ 

ISEQ subscript, first surface 
for this load-set 

Total number of surfaces calcu- 
lated thus far. 








Table 7. — ( Continued ) 


LAB 

DES 

ELED COMMON NAME: CVBMT2 

CRIPTION: VBMT subroutine linkage (temporary storage) 

NO. 

VARIABLE 

T 



DESCRIPTION 

1 

Q3 

R 

10000 


Temporary storage 
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Table 7. — (Continued) 


BLANK COMMON: 

NPLDS/PLDS Dynamic 

Storage 

DESCRIPTION: 

Variably dimensioned arrays 













NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

NM3 

I 

ISMAX 


Random Access Key for M 3 

2 

NM4 

I 

20*ISMAX 


Random Access Key for 

3 

NM5 

I 

20* ISMAX 


Random Access Key for m7 

4 

NC3 

I 

20* ISMAX 


Random Access Key for J 

5 

LIEOM 

R 

LBUF 


Buffer for IEOMLD 

6 

LIDISK 

R 

LBUF 


Buffer for IDISK 

7 

LNPTAP 

R 

LBUF 


Buffer for NPTAP 

8 

MERG 

R 

LBUF 


Buffer for MERGMB 

9 

LMERG 

R 

LBUF 


Random Access Index for MSRGMB 

10 

LMASS 

R 

LBUF 


Buffer for MASSTP 

11 

INODE 

I 

NODES 


Requested nodes 

12 

XMASS 

R 

NODES 

m . 
i 

Mass matrix 

13 

AREA 

R 

NODES 


Structural node areas 

14 

PMAT 

R 

NLOAD* 

NAERO 

p 

P - matrix 

15 

GEOM 

R 

3*NODES 

(X,Y,Z) s 

Structural Geometry 

16 

PZ 

R 

NODES* 

MXMODE 


<f> 2 of requested nodes 

17 

XYL 

R 

2*NAERO 

(x,y) l 

Local ( X , Y) coordinates 

18 

AREAL 

R 

NAERO 

[A] l 

Local aero panel areas 

19 

XKVAL 

R 

NK 

K 

Frequencies 

20 

FPL 

R 

NAERO* 

MXMODE 

1 1 

** 

r* 

Force coefficient matrix 
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Table 7. — (Continued) 


BLANK COMMON: 

VBMT Dynamic Storage 

DESCRIPTION: 

Variably 

dimensioned arrays 













NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

LJTAPE 

R 

LBUF 


FET buffer for JTAPE 

2 

LIEOM 

R 

LBUF 


FET buffer for IEOMLD 

3 

LIDISK 

R 

LBUF 


FET buffer for IDISK 

4 

LLTAP 

R 

LBUF 


FET buffer for LTAP 

5 

MERG 

R 

LBUF 


FET buffer for MERGMB 

6 

LINTMP 

R 

LBUF 


FET buffer for INTMP 

7 

ISEQ 

I 

ISMAX 


Sequential order of surface 
numbers 

8 

MVX 

I 

ISMAX 


Number of VX components for 
surface IS 

9 

MVY 

I 

ISMAX 


Number of VY components for 
surface IS 

10 

MVZ 

I 

ISMAX 


Number of VZ components for 
surface IS 

11 

MMXX 

I 

ISMAX 


Number of MXX components for 
surface IS 

12 

MMYY 

I 

ISMAX 


Number of MYY components for 
surface IS 

13 

MMZZ 

I 

ISMAX 


Number of MZZ components for 
surface IS 

14 

ISLOAD 

I 

ISMAX 


Sequential number of the first 
load for surface IS 

15 

LNODD 

I 

MXDUM 


Load number of defined dummy 
node (not the sequential 
number) 

16 

SCALD 

R 

MXDUM 

SCALED i 

Scale factor on dummy node 

17 

IDNOD 

I 

MXDUM 


Durrony nodes to sum to current 
surface 

18 

DAX 

R 

MXDUM 

bs d 

X coordinate of dummy node 

19 

DAY 

R 

MXDUM 

bbl d 

Y coordinate of dummy node 

20 

DAZ 

R 

MXDUM 

wl d 

Z coordinate of dummy node 
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Table 7. — (Continued) 


BLANK COMMON: 
DESCRIPTION:. 

VBMT Dynamic Storage (continued) 

Variably dimensioned arrays 









NO. 



VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

21 

Z 

R 

IZ 


Z override coordinates 

22 

NODEZ 

I 

IZ 


Node number for Z override 






coordinate 

23 

LMERG 

I 

(MXDUM+ 


MERGMB INDEX 




raLQAD) * 






(1+3NK) 



24 

NM3 

I 

MXLOAD 


Random access key number for M$ 

25 

NM4 

I 

MK*MXLOA 

3 

Random access key number for Mi 

26 

NM5 

I 

SJK*MXLOA 

3 

Random access key number for Ms 

27 

NC3 

I 

HK*MXLOA 

3 

Random access key number for $ 

28 

NMD 

I 

MXDUM* 


Random access key number for 




( 1+3NK) 


dummy nodes 

29 

XKVAL 

R 

NK 

K 

Frequencies 

30 

NROW 

I 

MXLOAD 


Number of rows in the matrix 






for this LOAD 

31 

THETX 

R 

MXLOAD 

e 

X 

THETAX associated with IL 

32 

THETY 

R 

MXLOAD 

e 

y 

THETAY associated with IL 

33 

THETZ 

R 

MXLOAD 

J 

e 

2 

THETAZ associated with IL 

34 

SAX 

R 

MXLOAD 

bs a 

X coordinate associated with 





1 

LNODD 

35 

SAY 

R 

MXLOAD 

BBL 

Y coordinate associated with 






LNODD 

36 

SAZ 

R 

MXLOAD 

WIj a 

Z coordinate associated with 






LNODD 

37 

INODS 

I 

NODES 

i 

Structural nodes for this 






surface 

38 

SCALS 

R 

NODES 

SCALES 

Scale factor for structural 






nodes 

39 

GEOM 

R 

6 * NODE? 

l 

x,y,z,6 ,9 ,0 of structural nodes 






2 X V z 

40 

INODA 

I 

NAERO 


Aero nodes for this surface 




Table 7. — ( Continued } 


BLANK COMMON : VBMT Dynamic Storage (continued) 

DESCRIPTION: Variably dimensioned arrays 


NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

description 

41 

SC ALA 

R 

NAERO 

SCALEA 

Scale factor for aero nodes 

42 

AXYZ 

R 

3 *NAERO 

BS.,BBL. , 

1 i' 

"Li 

X,Y,Z in the inertia axis for 
aero nodes (transformed from 
local X,Y coordinates) 

43 

XMA 

R 

NODES 
(or zero 

m . 

> 1 

M from the J-MATRIX 

44 

XMX 

R 

NODES 
(or zero 

( m i x i 

MX from the J-MATRIX 

45 

XMY 

R 

NODES 
(or zero 

m.y . 

MY from the J-MATRIX 

46 

XMZ 

R 

NODES 
(or zero) 

m . z . 
i i 

MZ from the J-MATRIX 

47 

XIXX 

R 

NODES 
(or zero) 

I 

XXI 

I XX from the J-MATRIX 

48 

XI YY 

R 

NODES 
(or zero) 

I 

yyi 

IYY from the J-MATRIX 

49 

XIZZ 

R 

NODES 
(or zero) 

I 

zzi 

IZZ from the J-MATRIX 

50 

XMXY 

R 

NODES 
(or zero) 

m . x . y . 
l i i 

MXY from the J-MATRIX 

51 

XMXZ 

R 

NODES 
(or zero) 

m . x . z . 

ill 

MXZ from the J-MATRIX 

52 

XMYZ 

R 

NODES 
(or zero) 

m . y . z . 
i*i i 

MYZ from the J-MATRIX 

53 

PX 

R 

STODES * 
4XMODE 

♦ 

X 

d> from INTERP 

y x 

54 

PY 

R 

tfODES * 
MXMODE 

>1 

<p^ from INTERP 

55 

PZ 

R 

NODES* 

MXMODE 

♦ 

z 

cf> from INTERP 
z 

56 

PPX 

R 

NODES* 

MXMODE 

<t> 

ex 

i > 0 from INTERP 
X 


83 




Table 7 . — (Concluded) 


BLANK COMMON : VBMT Dynamic Storage 

DESCRIPTION: Variably dimensioned arrays 


NO. 

VARIABLE 

T 

DIM. 

57 

l 

I 

PPY 

R 

NODES* 

MXMODE 

58 

PPZ 

R 

NODES* 

MXMODE 

59 

FPL 

R 

2*NAERO* 

MXLOAD 


DESCRIPTION 


k 0 from INTERP 


from INTERP 


Force coefficient matrix from 
EOM 









4.0 EXTENT OF CHECKOUT 


Each module of L218 (AVD, NPLDS/PLDS, and VBMT) was checked out with 
preliminary standalone data and then final verification test data, exercising the 
options indicated in tables 8 through 10. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1977 
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Table 9. — NPLDS/PLDS Checkout Summary 






7 

X 

X 

X 

X 

X 
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